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Impact of three years of data from the Wilkinson Microwave Anisotropy Probe on 
cosmological models with dynamical dark energy 
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The first three years of observation of the Wilkinson Microwave Anisotropy Probe (WMAP) have 
provided the most precise data on the anisotropics of the cosmic microwave background (CMB) to 
date. We investigate the impact of these results and their combination with data from other astro- 
physical probes on cosmological models with a dynamical dark energy component. By considering a 
wide range of such models, we find that the constraints on dynamical dark energy are significantly 
improved compared to the first year data. 
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INTRODUCTION 

Observations of type la supernovae (SNe la) [J, i, 
structure formation (LSS) [SL |j| and the cosmic mi- 
crowave background (CMB) [H, 0] all agree on an ac- 
celerated expansion of our Universe. This rather unex- 
pected phenomenon can be explained by modifying 4-D 
gravity @, Q or adding a new component to the total 
energy momentum tensor. The simplest such component 
is a cosmological constant. It fits all current observations 
flawlessly and has a simple interpretation in terms of a 
vacuum energy. Yet, its observed value is 120 orders of 
magnitude off from the naive estimate A ~ Mp, where 
Afp is the reduced Plank mass. The coincidence between 
this minute dark energy contribution and the observed 
energy density of matter is rather puzzling. If not given 
by chance via some sort of anthropic principle, it neces- 
sitates a mechanism that explains this coincidence. An 
immediate possibility is a coupling [lol . 11 . 1^ between 
(dark) matter and dark energy (though there might be 
problems due to quantum effects (isjV 
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Another solution is an attractor behavior 
of dark energy that leads to an almost constant ratio be- 
tween the fractional energy density ^d{z) of dark energy 
and the species otherwise dominating the expansion, i.e. 
photons and neutrinos during radiation domination and 
matter during matter domination. Coincidentally, such 
an attractor behavior corresponds to a scalar field with 
exponential potential that arises in string theories and 
when solving the cosmological constant problem from 
the point of view of dilatation symmetry [l4|. The non- 
vanishing ^d{z) at higher redshifts alleviates the problem 
of explaining the coincidence of matter and dark energy 
today « rij. Instead of fine tuning A to many or- 
ders of magnitude, the tuning needed is of the order of 
10""^. However, the tuning needed for such early dark en- 
ergy cosmologies increases the less dark energy there is 
at earlier times. A detection of early dark energy, on the 
other hand, would give crucial hints to fundamental laws 
of nature. The aim of this study therefore is to investi- 
gate the implications of the three year data of WMAP 



on dynamical dark energy models in general and their 
respective fractions of early dark energy. 

In view of the theoretical uncertainties many differ- 
ent techniques have been employed in the analysis of the 
dark energy, ranging from atttempts to reconstruct the 
potential of a scalar field dark energy (e.g. Ref. 17[) 
to the principal component approach of Ref. [l3|- We 
consider the redshift dependence of the fractional dark 
energy ^d{z) as a free function to be "measured" by ob- 
servation. We investigate in this note various parameteri- 
zations and an interpolated model. The possible coupling 
between dark energy and dark matter is neglected in this 
study. 



OBSERVATIONAL TESTS 

Dark Energy influences the expansion history of our 
Universe. In particular, the age t^, conformal horizons 
of today tq and at last scattering r/s and the sound hori- 
zon Ts at last scattering are modified. The effects can 
be understood analytically 1^ 2^, 21 1 using an effective 
description in terms of weighted averages relevant for the 
epoch of last scattering 



and structure formation 



r!d(r) 
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rJsf = [In atr - In acq] ^/ fid (a) din a, (2) 



In atr 
In ac, 



where at,- ~ 1/3. The effect of dark energy on the CMB 
is twofold. Through the modified expansion history, it 
changes the acoustic scale I a- In addition, it leads to 
a decay of the gravitational potential that is seen as an 
integrated Sachs- Wolfe contribution and a suppression 
of fluctuations on small scales |24| . The suppression of 



growth can be understood by looking at the equation 
of motion for cold dark matter perturbations inside the 
horizon, where the dark energy fluctuations are negligible 
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Here the derivative is with respect to conformal time r. 
In a matter Universe, fl^ = 1 and the solution is 6m oc a. 
With dark energy present, < 1 and the growth of 
structure slows down according to the solution of ([3]) [2^ 



dm OC a 



[V25-240d-l]/4 ^ l-3nd/5 
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This suppression starts as soon as a mode enters the 
horizon. As 6m cannot grow during radiation domina- 
tion, this leads to a red tilt of the CMB and matter 
power spectra up to the scale of the mode entering just at 
matter-radiation equality fceq„. All modes with k > /ceg„ 
have been inside the horizon before equality and are sup- 
pressed by the same factor. Hence, early dark energy 
mimics to some extend a running spectral index with the 
important difference that the running stops at kequ- AH 
in all, the suppression leads to a smaller erg compared to 
a A-CDM universe according to [3| 
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where r is the conformal horizon today and ui is a suit- 
ably defined average equation of state of dark energy. As 
a rough rule of thumb, an increase of f2sf by 10% leads 
to a decrease of erg by 50%. In the following numerical 
analysis, we will not use constraints on the overall nor- 
malization of the power spectrum, i.e. we marginalize 
over the bias of 2dF and SDSS. As we will see, the data 
we use nevertheless constrains fist- 

In contrast to linear growth, non-linear structure for- 
mation is enhanced in early dark energy cosmologies. 
The density contrast 6c corresponding to a collapsed 
structure is lower than in A-CDM |21|]. As the abundance 
is exponentially sensitive to 6c, the cluster abundance is 
considerably higher for a given ag, as compared to A- 
CDM. In particular, the abundance of clusters at higher 
redshift drops more slowly than in A-CDM, which will 
soon be probed by gravitational lensing and Sunyaev- 
Zel'dovich surveys. 



INVESTIGATED MODELS 

Constraints on the values of cosmological parameters 
are always model dependent. For this analysis, we there- 
fore select dark energy models with a large variety of 
different and in part opposite physical properties. This 
approach allows to identify the model dependencies of 
the best fit ranges for the standard cosmological param- 
eters, and their sensitivity to a change of the underlying 
dark energy behavior as compared to the A-CDM model. 



The first of these dark energy models is a leaping kinetic 
term model ( "LKT" ) , where a change of the kinetic term 
of a scalar field at late times leads to acceleration 2^ . In 
addition, we consider two models described by a param- 
eterization of the evolution of the dark energy fraction 
rid (a)- In one of these models 25|, Od rather slowly re- 
laxes from today's f2j to an asymptotic early-time value 
of and Od can be written as Od = r2d(wo,fi*). The 
other parameterization [2^] , which is a function of wq and 
rijj exhibits a faster variation of fid and is characterized 
by a minimum amount of dark energy f2^ throughout all 
cosmological epochs. This essentially fixes fid to fl'^ from 
early times until redshifts of a few. We also include two 
models which parameterize the equation of state i/;(a) of 
dark e nerg y. The versatile parameterization of Bassett 
et al. [27[, generalized by Corasaniti and Copeland in 
Ref. [28| ("C&C"), has four parameters for the equation 
of state, namely its value today wq, its value during the 
matter dominated era Wm as well as the scale factor a™ 
and width A™ of the transition between these values, 
so that w — ^(wo, Wrn, a™, A„j). In addition, we con- 
sider the simple parameterization w{a) — wq + will — a) 
[2^,[33| frequently used in the literature. Finally, we ana- 
lyze a model with J7d linearly interpolated in In a between 
values at z = 1, 3, 10, 100 and z = 1100, leaving con- 
siderable freedom for the variation of lid (a) at the cost 
of a rather large number of model parameters. 

In addition to their respective dark energy parame- 
ters, all models depend on the standard cosmological 
parameters: the present matter energy fraction f2m and 
baryon energy fraction f2b, the Hubble parameter h, opti- 
cal depth r, scalar spectral index and the initial scalar 
amplitude As, which we took into account using the ob- 
servationally relevant combination ln{10^'^ As) — 2t. We 
chose flat priors on all parameters. 

The equation of state of the LKT and the r2d(wo, il*)- 
models was not allowed to cross the cosmological con- 
stant boundary oi w = —1, and their fluctuations were 
treated like scalar field perturbations. The other models 
were allowed to cross w = —1. For these models, the 
speed of sound = 6p/6p was fixed at = 1, so that 
the perturbation equations for the pressure fluctuations 
remain well-defined even at the crossing. For the interpo- 
lated and the Q{wo,^^) models, this procedure was only 
adopted when the equation of state was close to the cross- 
ing, and the fluctuations were treated as scalar field per- 
turbations everywhere else, so that these models had the 
usual scalar field perturbations during almost the entire 
evolution. This treatment is necessary because a sing le 
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scalar field cannot traverse w = —1 

These models were compared to two different sets of 
data. Set I are the WMAP 3-year d ata, and set II consists 



of WMAP i5i]^OOMERANG'03 [36| , VSA [37|, CBI 



i, ACBAR [H for the CMB plus 2DF ^ and SDSS 
3 for LSS and SNe la data [JU combined. We omitted 
the baryon acoustic oscillation data [3§| , as it is currently 
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FIG. 1: Monte-Carlo results for the cosmological parameters. 
The left-hand side shows the la confidence intervals for the 
comparison of the models to the 3-year data of WMAP, the 
right-hand side corresponds to the results for the analysis with 
the combined set II. The shaded regions depict the la inter- 
vals for the A-CDM analysis from WMAP for the WMAP 
data only and the "all" set, respectively. The crosses show 
the value of the respective parameter for the best fit model in 
the Monte-Carlo chains. 



not as sensitive as CMB and LSS in constraining early 
dark energy [i^ . 



RESULTS 

The results are summarized in Figure [U which yields 
both the constraints on each of the standard parame- 
ters for all six of the models as well as their scatter. Also 
shown are the values for Qgf and as- As a first result from 
the Monte-Carlo analysis, we find that cosmological mod- 
els with a dynamically evolving dark energy component 
fit the data as well as A-CDM, but not better. Secondly, 
the constraints on the standard cosmological parameters 
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FIG. 2: One-sigma confidence intervals for the fraction of dark 
energy at different redshifts for the interpolated model. Red is 
the result obtained by comparing to the WMAP 3-year data, 
blue (lower) corresponds to the combined set II as described 
in the text. The black line shows the evolution of Qa for a 
A-CDM universe with the WMAP 3-year best fit values for 
the standard cosmological parameters. 



are well compatible with the results found by WMAP 
for a A-CDM cosmology, almost irrespective of the as- 
sumed behavior of the dark energy. Models that allow 
for a significant amount of dark energy at early times do 
have, however, a few significant features. Most promi- 
nently, they have a lower erg from the linear analysis, as 
expected. In the light of the rather low scalar spectral in- 
dex Us found by WMAP, one might have suspected that 
dynamical dark energy models would allow for a scalar 
spectral index that would be closer to n^. = 1, which 
corresponds to a scale invariant spectrum of the initial 
fluctuations. However, all investigated dynamical dark 
energy models show a preference for tIs < 1, extending 
the result found by WMAP to a wider range of cosmo- 
logical models. 

The one sigma bounds on the dark energy fraction for 
the interpolated model are shown in Figure [5) It is ap- 
parent from this plot that, while the preferred amount 
of dark energy in this model includes the values typical 
for A-CDM cosmologies, it also allows for much larger 
fractions of dark energy for redshifts z = 3 and higher. 
The evolution of the dark energy density ~ i^^h^ can 
be substantial and is not required to be monotonic in 
this model. This can lead to pronounced ISW contribu- 
tions. With rid as a free parameter at different z, the 
model can to some extend "manufacture" the shape of 
the TT-power spectrum at scales larger than the first 
peak, leading to a good fit to the WMAP-3 data. This, 
however, comes at the cost of severely suppressing ag and 
is in conflict with the data of set II. Furthermore, with 
WMAP-3 alone, the Hubble parameter is considerably 
less constrained for this model. This is due to the fact 
that with Qdiz) for z > 1 given by the model's parame- 
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ters independent of h, the acoustic scale Ia depends only 
very weakly on h. The relative independence of in 
one redshift bin from its neighboring values also leads to 
the Monte-Carlo algorithm finding a comparatively high 
number of well- fitting models with rather high Clgf. The 
interpolated models consequently have the highest values 
of Qsf, with the upper 2a limit reaching Hgf ^ 11% for 
WMAP-3 alone and 6.5% from the combined set II. 

The 2(7 upper bounds of the other models for Qgf from 
WMAP-3 alone vary from a very low 0.5%, caused by the 
choice of parameterization for the w{a) = wq + wi{l — a) 
model, to ~ 5% for the il{wo , Q^)-mode\ and ^ 7% for 
the model with = fldiwo,^^*). This represents a de- 
crease of about one to two percent compared to the first 
year data. For the combined set II, the upper bound 
decreases for all these models to less than about three 
percent, with the ri(wo, f^d)-model yielding the highest 
value of Osf ^ 4%. We recall that already a few percent 
of Qgf can have important effects on the abundance of 
nonlinear structure at high redshift [2l] |. 



CONCLUSIONS 

The three year data of WMAP was used to estimate 
cosmological parameters for a wide range of dynamical 
dark energy models. We have shown that this new data 
in combination with large scale structure data constrain 
the average amount of dark energy during the time of 
structure formation to f2st ^ 4% for a fair sample of 
dark energy models from the literature. We have also 
constructed a parameterization of lid (a) which linearly 
interpolates between the dark energy fraction at several 
redshift bins. This allows for a considerably higher frac- 
tion of Clsi ■ The analysis also shows that the values of the 
standard cosmological parameters for the dynamical dark 
energy models are well compatible with the values found 
by WMAP for A-CDM. The effect of dynamical dark en- 
ergy on (Tg and on the formation of nonlinear structure 
offer promising routes for further constraints on the time 
evolution of dark energy or a possible falsification of the 
A-CDM model. 

We acknowledge the use of the Legacy Archive for Mi- 
crowave Background Data Analysis (LAMBDA). Sup- 
port for LAMBDA is provided by the NASA Office of 
Space Science. 



* Electronic address: M.Doran®thphys.uni-he idelberg.de| 

^ Electronic address: G.Robbers@thphys.uni-heidelberg.de 

* Electronic address: C.Wetterich@thphys.uni-heidelberg.de 
] P. Astier et al, arXiv:astro-ph/0510447 

I] A. G. Riess et al. [Supernova Search Team 
Collaboration], Astrophys. J. 607, 665 (2004) 
,arXiv:astro-ph/0402512. . 



[3: 
[4: 
[6: 

[7] 

[8] 
[9] 

[lo; 
[11 

[12 

[13 

[14 
[15 
[16 

[17 

[18 

[19 

[20 

[21 

[22; 

[23; 
[24; 
[25; 

[26; 

[27 

[28; 

[29 

[30 
[31 

[32; 

[33 
[34; 

[35] 



W. J. Percival et al. [The 2dFGRS Collaboration], 
Mon. Not. Roy. Astron. Soc. 327, 1297 (2001) 



arXiv:astro-ph/0105252 



M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 
69 (2004) 103501 arXiv:astro-ph/0310723 . 
D. N. Spergel et al., arXiv:astro-ph/0603449 

A. C. S . Readhead et al, Astrophys. J. 609 (2004) 498 
[arXiv:a stro-ph/0402359 . 

J. H. G oldstein et al, Astrophys. J. 599, 773 (2003) 
[arXiv :astro-ph/0212517 . 

G. R. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. 
B 485, 208 (2000) arXiv:hep-th/0005016 . 
J. D. Bekenstein, Phys. Rev. D 70, 083509 
(2004) [Erratu m-ibid. D 71, 069901 (2005)] 
[arXiv:astro-"ph70 403694 . 

C. Wetterich, Astron . Astrophys. 301, 321 (1995) 
[arXiv:hep-th /9408025|. 

L. Amendola, Phys. Rev. D 62, 043511 (2000) 
arXiv:astro-ph/9908023 . 

L. Amendola, M. Quartin, S. Tsujikawa and I. Waga, 
arXiv:astro-ph/0605488 ^ 

M. Doran and J. Jaeckel, Phys. Rev. D 66, 043519 (2002) 
arXiv:astro-ph/0203018 . 
C. Wetterich, Nucl. Phys. B 302, 668 (1988) 

B. Ratra and P. J. Peebles, Phys. Rev. D 37, 3406 (1988) 
R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. 
Lett. 80, 1582 (1998) 

M. Sahlen, A. R. Liddle and D. Parkinson, 
|arXiv:astro-ph/0610812| 

DT Huterer and G. Starkman, Phys. Rev. Lett. 90, 031301 
(2003) arXiv:astro-ph/0207517 . 

M. Doran, M. J. Lilley, J. Schwindt and C. Wetterich, 

Astrophys. J. 559, 501 (2001) ;arXiv:astro-ph /OO 1 2 1 39) . 

M. Doran, J. M. Schwindt and C. Wetterich, Phys. Rev. 

D 64, 123520 (2001) arXiv:astro-ph/0107525j. 

M. Bartelmann, M. Doran and C. Wetterich, 

arXiv:astro-ph/0507257 

R. R. CaldweU, M. Doran, C. M. Mueller, G. Schae- 
fer and C. Wetterich, Astrophys. J. 591 (2003) L75 
arXiv:astro-ph/ 0302505. . 

P. G. Ferreira an d M. Joyce, Phys. Rev. D 58, 023503 
(1998) arXiv:astro-ph/9711102 . 

A. Hebecker and C. Wetterich, Phys. Lett. B 497, 281 
(2001) arXiv:hep-ph/0008205 . 

C. Wetter ich, Phys. Lett. B 594, 17 (2004) 
[arXiv:astro -ph/0403289 . 

M. Dor an and G. Robbers, JCAP 0606, 026 (2006) 
[arXiv:a stro-ph/0601544 . 

ST Bassett, M. Kunz, J. Silk and C. Un- 
Mon. Not. Roy Astron. Soc. 336, 1217 (2002) 



B. 

garelli, 
[arXiv:a stro-ph/0203383 . 

P. S. Corasaniti and E. J. Copeland, Phys. Rev. D 67, 
063521 (2003) arXiv:astro-ph/0205544 . 
M. Chevallier and D. Polarski, Int. J. Mod. Phys. D 10, 
213 (2001) arXiv:gr-qc/0009008 . 
E. V. binder, Phys. Rev. Lett. 90, 091301 (2003). 
A. Vikman, Phys . Rev. D 71, 023515 (2005) 
[arXiv:astro-ph/04071 07 . 
W Hu, Phys. Rev. D 71 



(2005) 047301 



[arXiv:astro-ph /0 410680). 

G. Huey, arXiv:astro-ph/0411102 

R. R. Caldwell and M. Doran, Phys. Rev. D 72 (2005) 
043527 arXiv:astro-ph/0501104 . 

G. B. Zhao, J. Q. Xia, M. Li, B. Feng and X. Zhang, Phys. 



5 



Rev. D 72, 123515 (2005) "arXiviastro-ph /0507482] . 
[36] C. J. MacTavish et al, arXiv:astro-ph/0507 503 
[37] C. Dickins on et al, Mon. Not. Roy . Astron. Soc. 353 

(2004) 732 [arXiv:as tro-ph/0402498' . 
[38] C. 1. Kuo et al . [AC BAR collaboration], Astrophys. J. 

600, 32 (2004) !arXiv:astro-ph/0212289j . 



[39] D. J. Eisenstein et al, Astrophys. J. 633, 560 (2005) 
,arXiv:astro-ph/0501171_. 

and 



[40] M. Doran, S. S tern 
|arXiv:astro-pIi/060"9075| 



E. Thommes, 



